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Abstract — This work presents joint iterative power allocation 
and interference suppression algorithms for spread spectrum net- 
works which employ multiple hops and the amplify-and-forward 
cooperation strategy for both the uplink and the downlink. We 
propose a joint constrained optimization framework that con- 
siders the allocation of power levels across the relays subject to 
individual and global power constraints and the design of linear 
receivers for interference suppression. We derive constrained 
linear minimum mean-squared error (MMSE) expressions for 
the parameter vectors that determine the optimal power levels 
across the relays and the linear receivers. In order to solve 
the proposed optimization problems, we develop cost-effective 
algorithms for adaptive joint power allocation, and estimation 
of the parameters of the receiver and the channels. An analysis 
of the optimization problem is carried out and shows that the 
problem can have its convexity enforced by an appropriate choice 
of the power constraint parameter, which allows the algorithms 
to avoid problems with local minima. A study of the complexity 
and the requirements for feedback channels of the proposed 
algorithms is also included for completeness. Simulation results 
show that the proposed algorithms obtain significant gains in 
performance and capacity over existing non-cooperative and 
cooperative schemes. 

Index Terms — DS-CDMA networks, cooperative communica- 
tions, joint optimization, resource allocation, cross-layer design. 



I. Introduction 

Multiple collocated antennas enable the exploitation of the 
spatial diversity in wireless channels, mitigating the effects of 
fading and enhancing the performance of wireless communi- 
cations systems. Due to size and cost it is often impractical 
to equip mobile terminals or sensor nodes with multiple 
antennas. However, spatial diversity gains can be obtained 
when terminals with single antennas establish a distributed 
antenna array through cooperation |IT]-||3]. In a cooperative 
transmission system, terminals or users relay signals to each 
other in order to propagate redundant copies of the same 
signals to the destination user or terminal. To this end, the 
designer must employ a cooperation strategy such as amplify- 
and-forward (AF) p], decode-and-forward (DF) |]3], E] and 
compress-and-forward (CF) 

Prior work on cooperative multiuser direct-sequence code- 
division multiple-access (DS-CDMA) systems in interference 
channels has focused on problems that include the impact 
of multiple access interference (MAI) and intersymbol in- 
terference (ISI), the problem of partner selection H, ifTOl 
and the bit error rate (BER) llT2l . |fT3l . outage performance 
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analysis issues ifTTl . and the evaluation of the spectral effi- 
ciency II22I and the diversity gains p3| based on asymptotic 
results. The main motivation for cooperative relaying with 
DS-CDMA systems is to increase the capacity, reliability and 
the interference suppression capability of these networks ifTOl . 
ifTTl . II22I . II23I . Recent contributions in the area of cooperative 
communications have considered the problem of resource 
allocation ID, Q in multi-hop time-division multiple access 
(TDMA) systems and MIMO systems [8|, (H. Related work 
on DS-CDMA system has focused on adaptive modulation 
|fT4l|. power and rate allocation |[T8l . ||20l and scheduling 
f21|. In the literature, there has been no attempt to jointly 
consider the problem of power allocation and interference 
mitigation in cooperative multiuser DS-CDMA systems so 
far. This problem is of paramount importance in cooperative 
wireless ad-hoc and sensor networks llT4l - llT8l that utilize DS- 
CDMA systems. These networks require multiple hops to 
communicate with nodes that are far from the base station 
in order to increase their coverage [lT9l . Moreover, multi-hop 
cooperative relaying can substantially improve the interference 
suppression capabilities |[TOl . |[T2l . |[T3l . 

The goal of this paper is to devise a cross-layer optimization 
strategy to significantly increase the capacity, reliability and 
coverage of spread spectrum networks which employ multi- 
ple hops and the AF cooperation protocol. Specifically, the 
problem of joint resource allocation and linear interference 
suppression in multiuser DS-CDMA with a general number of 
hops is addressed. In order to facilitate the receiver design, we 
adopt linear multiuser receivers ll24l . Il25l which only require a 
training sequence and the timing. More sophisticated receiver 
techniques 1241 . Il26l . 1271 . l28l are also possible for situations 
with increased levels of interference. A joint constrained 
optimization framework that considers the allocation of power 
levels among the relays subject to individual and global power 
constraints and the design of linear receivers is presented. It 
should be noted that the proposed design with individual power 
constraints has been initially reported in ll29l . whereas the pro- 
posed design with both individual and global power constraints 
has been introduced in ||30l. Here, the proposed designs are 
described and investigated in further detail, more complete 
derivations along with analysis and simulations results are 
included. Linear MMSE expressions that jointly determine 
the optimal power levels across the relays and the linear 
receivers are derived. Adaptive least squares (LS) algorithms 
are also developed for efficiently solving the joint optimization 
problems and mitigating the effects of MAI and ISI, and 
allocating the power levels across the links. An analysis of 
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the optimization problem is conducted and shows that the 
problem can have its convexity enforced by an appropriate 
selection of the power constraint parameter This allows the 
algorithms to avoid problems with local minima. A study 
of the computational complexity and the requirements for 
feedback channels of the proposed algorithms is also included. 
The main contributions of this work can be summarized as: 

1) A joint constrained optimization framework for the alloca- 
tion of power levels among the relays subject to individual and 
global power constraints and the design of linear receivers; 

2) Constrained linear MMSE expressions for the power allo- 
cation and the design of linear receive filters; 

3) Recursive algorithms for estimating the channels, the power 
allocation and the receive filters; 

4) Convexity analysis of the proposed optimization problems; 

5) A study of the computational complexity and the require- 
ments for feedback channels of the proposed and existing 
algorithms. 

The rest of this paper is organized as follows. Section II 
describes a cooperative DS-CDMA system model with multi- 
ple relays. Section III is devoted to the problem formulation 
and the constrained linear MMSE design of the interference 
mitigation receiver and the power allocation. The proposed 
LS algorithms for the estimation of the receive filter, the 
power allocation and the channels subject to a global and 
individual power constraints are developed in Sections IV and 
V, respectively. Section VI is devoted to the analysis of the 
computational complexity and feedback requirements of the 
proposed algorithms. Section VI presents and discusses the 
simulations and Section VII draws the conclusions. 

II. Cooperative DS-CDMA System and Data 
Models 

Consider a synchronous DS-CDMA system communicating 
over multipath channels with QPSK modulation, K users, N 
chips per symbol and L (L < N) as the maximum number 
of propagation paths for each link. The synchronous DS- 
CDMA systems is considered for simplicity as it captures 
most of the effects of asynchronous systems with a low delay 
spread ||25l , ||26l . We consider both uplink and downlink 
transmissions. The network is equipped with an AF protocol 
that allows communication in multiple hops using jir relays in 
a repetitive fashion. Therefore, we have Up ~ rir + I phases 
of transmission or hops and only one transmitter (source or 
relay) is active per phase, which increases the delay but also 
improves the coverage. The throughput is affected by the fact 
that there is an extra time slot per phase of transmission, 
however, there are situations for which the performance gains 
can offset the extra time slots and the throughput can be 
improved. Other cooperation protocols such as DF can be 
employed without significant modifications, however, the AF 
has been adopted for simplicity and due to its lower complex- 
ity for implementation ||3|. We assume that the source node 
or terminal transmits data organized in packets comprising P 
symbols, where there is a preamble with training symbols 
followed by a part with data symbols. We also assume that 
the packet contains a sufficient number of training symbols in 
the preamble for parameter estimation and that the network 
can coordinate transmissions and cooperation. The relays 



and destination terminals are equipped with linear receivers, 
which are synchronized with their desired signals. Since the 
focus of this work is on the resource allocation and linear 
interference mitigation, we assume perfect synchronization, 
however, this assumption can be relaxed to account for more 
realistic synchronization effects in the system. The proposed 
algorithms for power allocation and interference mitigation are 
employed at the receivers. A feedback channel is required to 
convey the power allocation parameters, which should cope 
with the channel variations. 
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Fig. 1 . (a) Uplink and (b) downlink of tlie cooperative DS-CDMA system. 

The cooperative DS-CDMA system under consideration is 
depicted in Fig. [T] The data model is described for the uplink 
in what follows. However, it should be remarked that the 
downlink data models can be obtained as a particular case of 
the uplink one. The received signals are filtered by a matched 
filter, sampled at chip rate and organized into M x 1 vectors 
fsd, fsTi and r^id, which describe the signal received from 
the source to the destination, the source to the relays, and the 
relays to the destination, respectively. 



Tsd = X! "'sdDkhsd,kbk + Vsd + ^sd, 

k=l 
K 

"^^^3 = X! a't^^Dkhsr^M^k + 'Hsr, + ^sr, , (1) 
k=l 
K 

Vr^d = X! '^jdDkhr^d.kVk ' + "Hr^d + "r,-<i, 



k=l 



where AI = N + L — \. The quantity hk[mj\ corresponds to 



3 



the transmitted symbol of user k, whereas fe^^*^ represents the 
symbol fe^^'^ processed at the relay j using the AF protocol. 
The amplitudes of the source to the destination, the source 
to the relay j and the relay j to the destination links for 
user k are denoted by a^^, a^^. and aj!,^, respectively. The 
vectors Ugd, ngr^ and rir-d represent the noise at the receiver 
of the destination and the relays. The vectors t]^^, t]^^ and 
T]^.^ denote the intersymbol interference (ISI) arising from 
the source to destination, source to relay j and relay j to 
destination links, respectively. The M x L matrix Dk has the 
signature sequences of each user shifted down by one position 
at each column that form 



(«) Ti: 



scheme: Terminal k 
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(b) Reception scheme: 
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stands for the 

signature sequence of user fc, the L x 1 channel vectors from 
the source to the destination, the source to the relay, and the 
relay to the destination are hsd.k, h^-d.k, hr-s.k, respectively. 
By stacking the data vectors in ^ (including the links from 
the relays to the destination) into a (n,, + 1)A/ x 1 received 
vector at the destination we have 

^.,iDhhxd.kh. 



Fig. 2. (a) Transmission and (b) reception sclieme of tlie cooperative DS- 
CDMA system. 



III. Problem Statement and Proposed MMSE 
Design 

This section states the problem of joint power allocation 
and interference suppression for a cooperative DS-CDMA 
network. Specifically, constrained optimization problems are 
formulated in order to describe the joint power allocation 
and interference suppression problems subject to a global 
and individual power constraints. The proposed linear MMSE 
designs are aimed at the destination, which is responsible 
for jointly computing the receiver parameters and the power 
allocation that is sent via a feedback channel to the source. 
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(3) 



By using the stacked received data from the source and the 
relays for joint processing and using i to denote the desired 
symbol in the transmitted packet and its received and relayed 
copies, we can rewrite the data in a compact form given by 

K 

r\i] =^Cknk\i]Bk[i]ak[i\+v[i]+'n[i], (4) 

k=l 

where the [n,. + 1)AI x (n^ + 1)L matrix Ck contains shifted 
versions of Dk as shown by 

Dk 



rXs^A. MMSE Design with a Global Power Constraint 

^rid The linear MMSE design of the power allocation of the 
: link, across the source, relay and destination terminals and 
inte ference suppression filters is presented here using a global 
power constraint. Let us express the received vector in (|4} in 
a more convenient way for the proposed optimization. The 
{rir + 1)A/ X 1 received vector can be written as 



r[i] = CT'HT\i\BT[i]aT[i] + r]\i\ 
K{nr 



(6) 



where the (ur + 
Ct = [Ci C2 ... 
the K{nr + l)L 
contains the channel 

K{nr + 1) X K{nr 
diag(6i[mi] b["'[m2].. 



diagonal matrix 

[rrinj ■ . .6k[toi 



- l)L matrix 
the signatures, 
matrix Ht [i] 

the links, the 
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(5) 



has the channel 



The {rir + 1)L x {ur + 1) matrix "Hkii] 
gains of the links between the source and the destination, 
and the relays and the destination. The (n^ + 1) x (n^ + 1) 
diagonal matrix JBfc[i] = diag(6fe b^f}'^ . . .b'^."'^) contains the 
symbols transmitted from the source to the destination (bk) 
and the rir symbols transmitted from the relays to the desti- 



1)M X 

Ck] contains all 
X K{nr + 1) 
gains of all 
+ 1 

contains the symbols transmitted from all the 
sources to the destination and from all the relays 
to the destination on the main diagonal, and the 
K{nr + 1) X 1 power allocation vector aT[i] = 
[ala[mi] al^^[m2] . . . ai_^_^[m„J . . . a^[mi] a^rf[m2] . . . a^^^^[m 
contains the amplitudes of all the links. 

Consider a joint MMSE design of the receivers for the 
K users represented by a {rir + 1)M x K parameter matrix 
W[i] = . ■ • , if/f [*]] and for the computation of the 

K{nr + 1) X 1 optimal power allocation vector ay.optH- This 
problem can be cast as 



nation (bl"^ . 
vector ak[i] =■ [a' 



6)^" ) on the main diagonal, the (ti^ + 1) x 1 
d^ °f amphtudes, the 



sd ^nd - ■ ■ ^r-„, 

{nr + l)M xl vector r][i] with the ISI terms and (n,. + l)Mx 1 
vector n[i] with the noise. A schematic that summarizes the 
transmission and reception schemes is depicted in Fig. |2] 



[W^opt, iT^opt 

subject to a;^[i]aT[i] = Pt 



arg mm 



E[\\b[i] - W"[i\r[i\ 



where the K x 1 vector b[i] = [bi[i], ■ ■ 
the desired symbols of the K users. 



(7) 

, represents 
The Unear MMSE 
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expressions for the parameter matrix Wopt and the vector 
ciT.opt can be obtained by transforming the above constrained 
optimization problem into an unconstrained one with the 
method of Lagrange multipUers ||39l which leads to 

Ct = E[\\b[i] - W"[i]{eTnT\i]BT[i\aT[i]+ri[i\+n[i\)\\ 

(8) 

Fixing aT[i], taking the gradient terms of the Lagrangian and 
equating them to zero yields 

W^opt = R 'Pen, (9) 

where the covariance matrix of the re- 
ceived vector is i? = £'['r['t]r-'^[i]] ~ 
CT'HT\i]BT\i]aT\i]al^[i\Bf^\i]'HT\i]CT + cr^/ and 
Pen = E[r[i]b"[i\] = E[CT-HT[i]BT[i]aT[i]b"[i\] is the 
{rir + 1)M X K cross-correlation matrix. The matrices R 
and Pc-H depend on the power allocation vector aT[«]. The 
expression for is obtained by fixing taking the 

gradient terms of the Lagrangian and equating them to zero 
which leads to 

where the K{nr + 1) x K{nr + 1) covariance matrix Raj, = 
E[B^ [{[Ht [ACt W[i\W^[i\CTnT[i\BT[i\] and the vector 
= E[B^[{\'HT[i]"C"W[i]b[i]] is a K{nr + 1) x 1 
cross-correlation vector The Lagrange multiplier Ay in the 
expression above plays the role of a regularization term and 
has to be determined numerically due to the difficulty of 
evaluating its expression. The expressions in (|9]l and (fTOt 
depend on each other and require the estimation of the channel 
matrix Ht\i\- Thus, it is necessary to estimate the channel and 
to iterate (|9]l and (fTOl i with initial values to obtain a solution. 
In addition, the network has to convey all the information 
necessary to compute the global power allocation including 
the filter Wopt- The expressions in (|9]l and (fTOl i require 
matrix inversions with cubic complexity ( 0{{nr + 1)M)'^) 
and 0{{K{nr + 1))'^), should be iterated as they depend on 
each other and require channel estimation. 

B. MMSE Design with Individual Power Constraints 

Here, the joint design of a linear MMSE receiver and 
the calculation of the optimal power levels across the relays 
subject to individual power constraints is presented. Consider 
an MMSE approach for the design of the receive filter Wk[{\ 
and the power allocation vector ak[{\ for user k. This design 
problem is posed as 

[■ujA;,opt,afc,opt] = arg min E[\bk[i] - it>f [i]r[i]|2] 

subject to af[z]aA;[i] = Pa. fe, A: = 1, 2, ...,K. 

(11) 

The expressions for the parameter vectors Wk[i] and ak[i\ can 
be obtained by transforming the above constrained optimiza- 
tion problem into an unconstrained one with the method of 



Lagrange multipliers ||39l . which leads to 

K 

Ck = E[\bk[i]-w^\i]{Y,CiHi\i\Bi\i\ai[i]+'n\i\+n[^)\^] 
1=1 

+ Afc(af Hafc[z] -PA,fc). k = l,2, ...,K. 

(12) 

Fixing ak[i], taking the gradient terms of the Lagrangian and 
equating them to zero yields 

Wk,opt = R^'PcH, k = l, 2, . . . , K, (13) 

where R = Y.k=i Cfc?^feWBfc Wofc Waf W-Bf ^Cf + 
a^I is the covariance matrix and Pq-h ~ -^[^feM'^'I*]] ~ 
Cfc^[i]afe[i] is the cross-correlation vector The quantities R 
and Pen depend on ak[i\. By fixing tOfc[i], the expression for 
ak [i] is given by 

flfe.opt = (fia, + AfeJ)-Va,, /c= 1, 2, i^, (14) 

where R^, = Eti M?if MCf t(;,H«;f WCfe^fcHB^H 
is the {rir + 1) x (n^ + 1) covariance matrix 
and the (n,. + 1) x 1 cross-correlation vector is 
P„ = E[bk[t]B^[mk[A''Cfwk[2]]. The expressions 
in ( fTSl ) and ( fT4l i have to be iterated as they depend on each 
other and require the estimation of the channel matrices 
Hkii]- The expressions in ( fT3T l and (fl4] l also require matrix 
inversions with cubic complexity ( 0{{{nr + 1)M)^) and 
0{{nr + 1)"^). In what follows, we will develop adaptive 
algorithms for computing opt, iffc.opt and the channels 
Hkli] for /s = 1, . . . , A' in an alternating fashion. 

IV. Proposed Joint Estimation Algorithms with a 
Global Power Constraint 

Here we present adaptive joint estimation algorithms to 
determine the parameters of the linear receiver, the power 
allocation and the channel with a global power constraint. The 
proposed joint power allocation and interference suppression 
(JPAIS) algorithms with a global power constraint (GPC) are 
simply called JPAIS-GPC, are suitable for the uplink of DS- 
CDMA systems and rely on LS-based estimation algorithms. 
The proposed algorithms are based on the idea of alternating 
optimization BDl . Il32l . in which the recursions for computing 
the parameters of interest are employed in cycles of iterations 
and over the received symbols. Note that more advanced 
algorithms ||33l - ll38l could also be considered in this context. 

A. Receiver and Power Allocation Parameter Estimation Al- 
gorithms 

Let us now consider the following proposed least squares 
(LS) optimization problem 

i 

[#W,aTW] - arg min Y.^'-'\\b[l] - 
subject to [{\aT[i] ~ Pt, 

(15) 

where a is a forgetting factor The goal is to develop a 
cost-effective recursive solution to ( flSt . To this end, we will 
resort to the theory of adaptive algorithms [|39| and derive 
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a constrained joint iterative recursive least squares (RLS) 
algorithm. This algorithm will compute W[i] and aT[i] and 
will exchange information between the recursions. The part 
of the algorithm to compute Wli] uses = R[i] ~ 



and is given by 



1 + Q;-ir-f^[i]*[i]r[i] ' 



= a-^^[i - 1] - a-^k\i\r"\i\^[i - 1], 

where the a priori estimation error is given by 

^ b[i] - w" [i - \\r[i]. 



(16) 

(17) 
(18) 

(19) 



The derivation for the recursion that estimates the power 
allocation clt [i] presents a difficulty related to the enforcement 
of the constraint and how to incorporate it into an efficient LS 
algorithm. This is because the problem in (fT5i incorporates a 
Lagrange multiplier (At) to ensure the global power constraint 
and the resulting system of equations cannot be solved with 
the aid of the matrix inversion lemma ||39l due to its structure. 
Our approach is to obtain a recursive expression by relaxing 
the constraint, to solve the system of equations with Gaussian 
elimination or with an adaptive conjugate gradient algorithm 
BTI . and then ensure the constraint is incorporated via a 
subsequent normalization procedure. 

The power allocation vector aT[i] is estimated via the 
solution of the following equation 



(20) 



where the K{nr + 1) x A' (n,. + 1) input data correlation matrix 
is given by 



= Qfi„^[i - 1] + Ut\i]u^ W, 



(21) 



where the K{nr + 1) x K matrix UT[i] = 
Bf^miij. [l]CTW[i\ and the K{nr + 1) x 1 cross-correlation 
vector is given by 



Put m = Paj 



i] + UT[m]- 



(22) 



The system of linear equations in (|20] | can be solved via 
Gaussian elimination with cubic complexity or more efficiently 
with quadratic complexity via an adaptive conjugate gradient 
algorithm BTI as described below 



W Par m - -""T m O'T [l\ , Cl[l\^V [l\ , 

axii + 1] = otH + "W-RdT 

V[i+ 1] = V[i] - Q![i]fiaTM'iHi 

+ 1] = v"[i + l]v[i + l]/{v"[i]v[t]), 

d[i + l] = v[i + 1] + I3[i + l]d\i], 

In order to ensure the global power constraint d^[i]dT[i] = 
Pt, we apply the following rule 

' , -1 



(29) 



(23) 
(24) 
(25) 
(26) 
(27) 
(28) 



drii] ^ \/Pt dT[i]{y a"[i]dT[i]) 



The algorithms for recursive computation of W[i] and aT[i] 
require estimates of the channel vector Ht[i\, which will 
also be developed in what follows. The complexity of the 
proposed algorithm is 0{{{nr + 1)A/)^) for calculating W^[z] 
and 0{{K{nr + 1))^) for obtaining otW- 

B. Channel Estimation with a Global Power Constraint 

We present a channel estimator that considers jointly all the 
K users and exploits the knowledge of the receive filter matrix 
W[i] and the power allocation vector dT[«]. Let us consider 
the received vector in (|6) and develop a channel estimation 
algorithm for Ht[i\- The proposed channel estimator can be 
derived from the following optimization problem 

i 

iiT[i] = arg min Va'-'||r[;] - CTHT[i\BT[i]aT[i\\\'^ . 

(30) 

The solution to the above optimization problem is given by 



(31) 



where the matrix inversion (C^Ct) ^ can be pre-computed 
and stored at the terminal of interest, and the /\ + 1)L x 
{rir + 1)L correlation matrix P<n^[i] is computed by the 
formula 



1=1 



(32) 



^PnAi-n+Cf^r[i]u^^[il 

where the K{nr + 1)L x 1 input data vector u-u^[i] for this 
recursion is 

u-Hrii] = BT[i\aT[i\, (33) 

and the inverse of the K{nr + l)L x K(nr + l)L matrix 
-R«T W = 121=1 '■'^'Ht W'^ILJ,^] is computed with the aid 



of the matrix inversion lemma fl39l as follows 



a 



-1 



1 + a [i]R-HT [« - ^u-Ht W 



(34) 

This algorithm jointly estimates the coefficients of the chan- 
nels across all the links and for all users subject to a global 
power constraint. Therefore, the complexity of the proposed 
RLS channel estimation algorithm is 0{{K{nr + 1)L)^). A 
summary of the main steps of the GPAIS-GPC algorithms is 
provided in Table I] 

TABLE I 

Summary of the JPAIS-GPC. 



1. Initialise parameters: W^[0], ot^O] and A. 
for i = 1 , . . . , P do 

2. Compute the receive filter W[i] using (18). 

3. Calculate the power allocation vector ax [i] using (25). 

4. Normalise aT[i] with (29) . 

5. Compute the channel estimate Jfy [i] using )3U . 
end for 

6. Transmit power allocation dy[i] to users. 
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V. Proposed Joint Estimation Algorithms with 
Individual Power Constraints 

In this section, we present adaptive joint estimation algo- 
rithms to determine the parameters of the linear receiver, the 
power allocation and the channel subject to individual power 
constraints. The proposed JPAIS algorithms with individual 
power constraints (IPC) are simply called JPAIS-IPC and, 
unlike the algorithms presented in the previous section, are 
more appropriate for the downlink and for distributed resource 
allocation, detection and estimation. Specifically, the JPAIS- 
IPC algorithms can be employed in a distributed fashion in 
which the linear receivers compute the individual receive filter, 
the power allocation and channel parameters, and then convey 
the power allocation parameters to the base station. 

A. Receiver and Power Allocation Parameter Estimation 

We develop a recursive solution to the expressions in (IT?t 
and ( fT4b using time averages instead of the expected value. 
The proposed RLS algorithms will compute Wk[i] and dfe[i] 
for each user k and will exchange information between the 
recursions. We fix ak[i] and compute the inverse of iR[i] using 
the matrix inversion lemma fi9\ to obtain ii»fc[i]. Defining 
= R[i] then we can obtain the recursions 



k[{\ = 



(35) 



^[i] = - 1] - a-^k[{\r"[i]^[i - 1], (36) 

Wk\i]=Wk[i-l] + mCk\i]. (37) 
where the a priori estimation error is 

S,k\i]=bk\i]^w^[i-l]r\i]. (38) 

The derivation for the recursion that estimates the power 
allocation requires a strategy to enforce the constraint on 
the individual power of user k, which circumvents the need 
to compute the Lagrange multiplier Afe. In order to develop 
the recursions for dfe[z], we need to compute the inverse 

of = ELi^f mCmcfri;fem«,f wc^fems.H. 

To this end, let us first define 4>aj. = Ra.^ [«] and proceed as 
follows: 



B. Channel Estimation with Individual Power Constraints 

We propose here an algorithm that estimates the channels 
for each user k across all links subject to individual power 
constraints and exploits the knowledge of the receiver filter 
Wkli] and the power allocation vector ak[i\. Let us consider 
the received vector in for the channel estimation procedure. 
A channel estimation algorithm can be developed to solve the 
following optimization problem 



for A: = 1, 2 



arff mm 



1=1 



'~'^^r[l]^ek-Hk[i]Bk[l]ak[l]\\\ 

(44) 



The solution to the above optimization problem is 



?i,W = (CfCfc) ^P^.WR. [z], for k^l,2, 



,K, 
(45) 

where the matrix inversion {C^Ck) can be pre-computed 
and stored at the base station, relays and mobile terminals, 
and the (n,. + 1)L x (n,. + 1)L correlation matrix -Ph^W is 
computed by the formula 



1=1 



(46) 



P^,[z-1]+Cfr[z]<[z], 



where the [rir + l)i x 1 input data vector u^^. [i] for this 
recursion is 

UhM"^ Bk[i\ak[i\ (47) 

and the inverse of the (n^ + \)L x {rir + 1)L matrix Rfi^ [i] 
is computed as follows 

(48) 

This algorithm estimates the coefficients of the channels of 
each user k across all the links subject to individual power 
constraints. The complexity is 0(((n,, + 1)L)'^). A summary 
of the main steps of the GPAIS-IPC algorithms is provided in 
Table I] 



(39) 
(40) 

ak[i] = ak[i] + k^MaM (41) 
where the a priori estimation error for the above recursion is 

Ui^ = - a^m^mfm^Wki^], (42) 

In order to ensure the individual power constraint 
aj^[i]ak[i] ~ PA,k, we apply the rule 



11, 



TABLE II 

Summary of the JPAIS-IPC. 



flfeW ^ a^[i]{^Ja^[i]a^[i]) . (43) 

The algorithms for recursive computation of Wk[i] and dfe[i] 
require estimates of the channel vector Hkii], which will 
also be developed in what follows. The complexity of the 
proposed algorithm is 0{{{nr + 1)M)^) for calculating Wk[i] 
and 0{{nr + 1)^) for obtaining dk[i]- 



1. Initialise parameters: lifclO], afc[0] and A. 
for each user k do 

for i = 1 , . . . , P do 

2. Compute the receive filter liifefi] using )37) . 

3. Calculate the power al loca tion vector ai^[i] using )4U . 

4. Normahse dfc[j] with j43j. 

5. Compute the channel estimate Hk[i] using 445t . 
end for with variable i. 

6. Transmit power allocation dj; [i] to base station. 



VI. Analysis and Requirements of the Algorithms 

In this section, we analyze the optimization problems de- 
scribed in Section III, and assess the requirements of the 
proposed JPAIS -GPC and JPAIS-IPC algorithms in terms of 
computational complexity and number of feedback bits. 
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A. Analysis of the Optimization Problems 

In this part, we develop an analysis of the optimization 
problems introduced in Section III, which can be extended to 
the algorithms introduced in Sections IV and V. Our approach 
is based on an algebraic manipulation of the original problem 
and the incorporation of the power constraints to illustrate 
the properties of the problem. The optimization problems 
considered in Section III constitute non-convex problems, 
however, it turns out that it is possible to devise a strategy 
that modifies the MSE-based optimization and enforces the 
convexity. This is corroborated by our studies that verify that 
the algorithms converge to the same solutions regardless of 
the initialization when this strategy is adopted. 

Let us consider the proposed optimization method in (|7) 
and examine the cost function 



C[W^],ak[i]) = E[\\b[i] -W'- \i\r\i\ 

K 



(49) 



k=l 



which must be minimized subject to [i]aT[«] — Pt and 
where [W^]fe[i] denotes the /c-column of the matrix W[i]. Let 
us now express x[i] = W^[i]r[i] as 



W"\i]r\i] 



W"[i]{CTnT[i]BT[i] arii] + W"[i]iv[i] + n[i]) 

m 



K 



k=l 



(50) 



where vr^ is a A' x 1 vector that contains a 1 in the fc-th 
position and zeros elsewhere. Now let us consider the joint 
optimization problem via a single parameter vector defined as 



(51) 



Using the above definition, we can rewrite the data symbol 
expression for user k as follows 



Xk [i\ 



K 

E 

fe=i 



K 

fe=l 



(52) 



where the (n.^ + + A') x {nr + l){M- 
the signal and noise components are 



A') matrices with 



U,\i] = 







(n^ + l)Kx(n,^ + l)K 

t\i] 



(nr + l)Mx(nr + l)M 
\nr + l)Kx(nr + l)K 



and 



(n,, + l)Afx(rv + l)(-fS'+A/)-l 
0(rir + l)/<'xl ^{nr + l)Ky.(nr + mK+M)-l 



If we now rewrite the cost function in (149) as a function of 
qk[i\, we obtain the equivalent function 



K 



fc=l 



(53) 



which is also subject to the global power constraint 
a;^[i]aT[«] = Pt- In order to evaluate the convexity of the 
optimization problem, we can verify if the Hessian {Hk) 
with respect to user k of C(qj.[i]) is positive semi-definite if 
m^Hkm > for all nonzero m e c'^^^^+i)('^i+K)xi gj^ 
Pol . Computing the Hessian HDl of the above cost function 
for the fcth user we obtain 



Hi 



' dq^m dq^m 

--E[{qfmT[^q,[^~bmUT[^ 
+ (qfWC/TWq,W-6feW)J7?W 



(54) 



By examining Hk, we verify that the second and fourth terms 
are positive semi-definite, whereas the first and the third terms 
are indefinite. At this stage, it is of central importance in 
the analysis to employ the constraint a^[i]ax[i] = Pt and 
perform some algebraic manipulations in the first and the third 
terms, which yield 



bmuTm^ 

(55) 



where /9y[i] = (aTa^[«])^aT, (O^ represents the pseudo- 
inverse and the operator 3fie{-} selects the real part of the 
argument. The above development shows that the optimization 
problem can have its convexity enforced by adjusting the 
power constraint Pt so that the following condition holds 



Pt > 



m"^e{E[{[W]^[t]nii]f3T[i])UT[i]] }m ' 



(56) 



The previous analysis can be conducted for the proposed 
JPAIS-IPC algorithm and leads to the following condition 

m"^e{E[{bl[i] - w^[2]j[i])UT[j]] }m 

TnH^e{E[{w^[i\Rk\ihk\i])UT\i]] }m ' 

where Uk[i] = Cknk[i]Bk[i], Uk+j[i] = Uk\i] 



Pk > 



(57) 



U,[2] = 



UM = 



'J(n,. + l)x(r: 

?fik[i] 

M 

0(n,. + l)xl 















(n,. + l)x(rv + l) 







1)M X 1 vector = Y.f^k Cini[i]Bi[i] + ri[i] + n[i], and 
7j,[i] = [aka^[i])'^ ak- These conditions have been verified 
by numerical experiments, which corroborate the analytical 
results. Therefore, the optimization problem can have its 
convexity enforced by an appropriate selection of the global 
and individual power constraints Pt and Pk, respectively. 

B. Computational Complexity Requirements 

We discuss here the computational complexity of the pro- 
posed and existing algorithms. Specifically, we will detail 
the required number of complex additions and multiplications 
of the proposed JPAIS-GPC and JPAIS-IPC algorithms, and 
compare them with interference suppression schemes without 
cooperation (NCIS) and with cooperation (CIS) ||T0|, HI] 
using an equal power allocation across the relays. Both uplink 
and downlink scenarios are considered in the analysis. In Table 
III we show the computational complexity required by each 



and 



the (n,. 
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table iii table v 

Computational complexity of algorithms with a global Computational complexity of the proposed JPAIS and 

POWER constraint . EXISTING ALGORITHMS. 





Number of operations per symbol 


Algorithm 


Recursions 


Parameter 


Additions 


Multiplications 








2(K + 1)M)2 


3(K + ^)Mf 


JPAIS-GPC (Uplink) 


W[i\, aT\i],HT[i\ 


W[i\ 


+2K{{nr + 1)M) 
-{Ur + l)M + 1 


+2K{{nr + l)M) 
+3(nr + 1)M + 1 


JPAIS-IPC (Downlink) 


Wkli], ak[i], Hk[i] 



3K{K{nr + 1)) 
+K{nr + 1){L-1) 
aT[i\ +K{M{nr + l)) 

+K{K{nr + 1)) 
+Q{K{nr + l)f 
+?,K{nr + 1) + n,. + 2 



H 



K{K{nr + l)) 
+'i{K{nr + l)f 
+ {K + L){K{nr + l)f 
-{K{nr + l)f 
+K{M{nr + l)L 



5(K(nr + l)Lf 
+5K{nr + 1)L 
+3 



+6K{nr + 1)L 
+ 1 



CIS (Uplink) 



aT[i] is fixed 



CIS (Downlink) 


Wk[l 


, ak[i] is 


fixed 


NCIS (Uplink) 


W[i\ 


with rir 


= 


NCIS (Downlink) 


Wk[i' 


with rir 


= 



TABLE IV 

Computational complexity of algorithms with 
individual power constraints. 



Number ot operations per symbol 



Parameter 


Additions 


Multiplications 


Wk[i] 


2{{nr + l)Mf 
+{nr + l)M 
+1 


3((n. + l)Mf 
+5(nr + 1)M 
+1 


ak[i] 


2{nr + lf 
+3{nr + 1) 
+M{nr + 1)L 
+ {nr + 1)L 
-3 


3(n, + if 
+7{nr + 1) 
+A/(n,. + 1)1/ 
+ {nr + 1)L 
+3 


Hk[i] 


2(K + l)i)2 
+5M(n,. + 1)L 
-5(nr + 1) + 3 


6(K + 1)L)2 
+A/(n^ + l)L 
+4:(nr + 1) + 1 



complexity is detailed in Tables III and IV. 

Complexity of RLS Algorithms 



i 10' 



10' 




-I- JPAIS-GPC ; 
-* — JPAIS-IPC 

B— CIS 

NCIS 

35 



recursion associated with a parameter vector/matrix for the 
JPAIS-GPC, which is more suitable for the upUnk. 

In Table IV we describe the computational complexity 
required by each recursion associated with a parameter vector 
for the JPAIS-IPC algorithm, which is suitable for both the 
uplink and the downlink. A noticeable difference between the 
JPAIS-GPC and the JPAIS-IPC is that the latter is employed 
for each user, whereas the former is used for all the K users 
in the system. Since the computation of the inverse of R[i] is 
common to all users for the uplink in our system, the JPAIS- 
GPC is more efficient than the JPAIS-IPC when the latter is 
computed for all the K users. 

The recursions employed for the proposed JPAIS-GPC and 
the JPAIS-IPC are general and parts of them are used in 
the existing CIS and NIS algorithms. Therefore, we can use 
them to describe the required computational complexity of 
the existing algorithms. In Table V we show the required 
recursions for the proposed and existing algorithms, whose 



Fig. 3. Computational complexity in terms of the number of complex 
multiplications of the proposed and existing schemes for the uplink. 

In Fig. [3] we illustrate the required computational com- 
plexity for the proposed and existing schemes for different 
number of relays (n^). The curves show that the proposed 
JPAIS-GPC and JPAIS-IPC are more complex than the CIS 
scheme and the NCIS. This is due to the fact that the power 
allocation and channel estimation recursions are employed. 
However, we will show in the next section that this additional 
required complexity (which is modest and equivalent to an 
additional 20 — 25% cost as compared to the CIS scheme) can 
significantly improve the performance of the system. 

C. Feedback Channel Requirements 

The proposed JPAIS algorithms require feedback signalling 
in order to allocate the power levels across the relays. In 
order to illustrate how these requirements are addressed, we 
can refer to Fig. |4] which depicts the structure for both the 
data and feedback packets. The data packet comprises a 
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preamble with a number of training symbols (A^tr), which 
are used for parameter estimation and synchronization, and 
the transmitted data symbols (A'data)- The feedback packet 
requires the transmission of the power allocation vector 
for the case of the JPAIS-GPC algorithm, whereas it requires 
the transmission of for each user for JPAIS-IPC. A typical 
number of bits ni, required to quantize each coefficient of the 
vectors ar and via scalar quantization is rif, = 4 bits. More 
efficient schemes employing vector quantization B2l . ||43l and 
that take into account correlations between the coefficients are 
also possible. 

For the uplink (or multiple-access channel), the base station 
(or access point) needs to feedback the power levels across 
the links to the K destination users in the system. With the 
JPAIS-GPC algorithm, the parameter vector with (n^ + 
l)Kni, bits/packet must be broadcasted to the K users. For the 
JPAIS-IPC algorithm, a parameter vector with (n^ + l)nf, 
bits/packet must be broadcasted to each user in the system. In 
terms of feedback, the JPAIS-IPC algorithm is more flexible 
and may require less feedback bits if there is no need for a 
constant update of the power levels for all K users. 



Data Packet Structure 



Feedback Packet Structure 
(JPAIS-GPC) 



Feedback Packet Structiue 
(JPAIS-IPC) 



{n,. + l)Knj, bits/feedback packet 



K times [n,- + l)Kni, bits/feedback packet 



Fig. 4. Proposed structure of the data and feedback packets. 

For the downlink (or broadcast channel), the K users must 
feedback the power levels across the links to the base station. 
With the JPAIS-GPC algorithm, the parameter vector ax with 
(rir + l)Knb bits/packet must be computed by each user and 
transmitted to the base station, which uses the ax vector 
coming from each user An algorithm for data fusion or a 
simple averaging procedure can be used. For the JPAIS-IPC 
algorithm, a parameter vector with (n^ + l)rib bits/packet 
must be transmitted from each user to the base station. In terms 
of feedback, the JPAIS-IPC algorithm requires significantly 
less feedback bits than the JPAIS-GPC in this scenario. 

VII. Simulations 

In this section, a simulation study of the proposed JPAIS 
and existing algorithms is carried out. The first existing scheme 
that is considered in the comparisons is a linear interference 
suppression technique that only takes into account the source 
to destination links and does not consider the contribution of 
the relays. This scheme is denoted non-cooperative interfer- 
ence suppression (NCIS) and corresponds to a Unear receive 
filter designed according to the MMSE criterion or computed 
with an RLS algorithm for each user The second existing 
scheme is denoted cooperative interference suppression (CIS) 
and processes the signals arriving from the source and the 
relays using an equal power allocation across the relays for 
each user The CIS scheme employs a linear receive filter 



designed according to the MMSE criterion or adjusted with 
an RLS algorithm, and the entries of the power allocation 
parameter vectors are equal (equal power allocation). We 
first evaluate the bit error ratio (BER) performance of the 
proposed JPAIS-GPC and JPAIS-IPC algorithms and compare 
them with the NCIS and the CIS schemes. We consider a 
DS-CDMA system with randomly generated spreading codes 
with a processing gain N = 16. The noise samples at the 
receivers of the relays and the destination are drawn from 
zero mean complex Gaussian random variables with variance 
a^. The fading channels (that can be time- varying or time- 
invariant) are generated using a random power delay profile 
with gains taken from a complex Gaussian variable with unit 
variance and mean zero, i = 3 paths spaced by one chip, and 
are normalized for unit power. The time-varying channels are 
generated according to Clarke's model B4l . which is param- 
eterized by the normalized Doppler frequency fdT, where fd 
is the Doppler frequency and T is the inverse of the symbol 
rate. The power constraint parameter Pam is set for each 
user so that one can control the SNR (SNR = Pam/'^'^) and 
Pt = KPa.u, whereas the power distribution of the interferers 
follows a log-normal distribution with associated standard 
deviation of 3 dB. We adopt the AF cooperative strategy 
with repetitions and all the relays and the destination terminal 
are equipped with linear MMSE or adaptive receivers. Note 
that the noise amplification of the AF protocol is considered 
O. The receivers have either full knowledge of the channel 
and the noise variance ( MMSE design) or are adaptive and 
estimate all the required coefficients and the channels using the 
proposed and existing algorithms with optimized parameters. 
For the JPAIS algorithms employing the MMSE expressions 
we employ 2 iterations per packet (when the channels are 
time-invariant) or per symbol (when the channels are time- 
varying) for the design of the parameter vectors, whereas for 
the adaptive versions we use only 1 iteration per update. The 
JPAIS algorithms are used at the destination and employ a 
feedback channel to send the power allocation vector to the 
source, whereas the relays are equipped with conventional 
linear MMSE or adaptive receivers. We employ packets with 
1500 QPSK symbols and average the curves over 1000 runs. 
For the adaptive receivers, we provide training sequences with 
iVtr = 200 symbols placed at the preamble of the packets. 
After the training sequence, the adaptive receivers are switched 
to decision-directed mode. 

We first consider the proposed JPAIS method with the 
MMSE expressions of (|9]l and ( fTOl i using a global power con- 
straint (JPAIS-GPC), and ([T3]l and (O with individual power 
constraints (JPAIS-IPC). We compare the proposed scheme 
with a non-cooperative approach (NCIS) and a cooperative 
scheme with equal power allocation (CIS) across the relays 
for Tir = 1,2 relays. The results shown in Fig. |5] illustrate the 
performance improvement achieved by the proposed JPAIS 
scheme and algorithms, which significantly outperform the 
CIS and the NCIS techniques. As the number of relays is 
increased so is the performance, reflecting the exploitation 
of the spatial diversity. In the scenario studied, the proposed 
JPAIS-IPC approach can accommodate up to 3 more users 
as compared to the CIS scheme and double the capacity 
as compared with the NCIS for the same performance. The 
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N=16, K=6 users, P-1500 



N-16, SNR= 12 dB, P-1500 




N=16, SNR - 12 dB, K=6 users, P-1500 



— e NCIS(n -2) 

—s — CIS(n -2) 

— K JPAIS-IPC(n.2) 

JPAIS-GPC(n.2) 




Number of users 



500 1000 
Number of received symbols 



Fig. 5. BER versus SNR and K for the optimal linear MMSE detectors. 
Parameters: \t = Aj. = 0.025. 



proposed JPAIS-GPC is superior to the JPAIS-IPC and can 
accommodate up to 2 more users than the JPAIS-IPC, while 
its complexity is higher. Equivalently, the BER versus SNR 
curves show that the proposed JPAIS scheme and algorithms 
can obtain a higher diversity order than the existing schemes, 
and save up to 2 dB in SNR for the same BER as compared 
with the existing techniques. The results in Fig. |5] suggest 
that the JPAIS-GPC algorithms are more suitable than the 
JPAIS-IPC algorithms for the uplink and situations with a 
high SNR and a large number of users. For the downlink and 
situations with low SNR and a small number of users, the gains 
of the JPAIS-GPC algorithms over the JPAIS-IPC algorithms 
are more very significant, suggesting that the latter are more 
suitable in these scenarios. The reason for the improved 
performance of the JPAIS algorithms over the existing schemes 
is that they jointly optimize the linear receive filter parameters 
and the power allocation, better exploiting the degrees of 
freedom at both the transmitter via power allocation and at 
the receiver with linear interference suppression to mitigate the 
interference. This approach allows a more effective reduction 
of the MSE and an improvement in the BER. 

The second experiment depicted in Fig. |6] shows the BER 
performance of the proposed adaptive algorithms (JPAIS) 
against the existing NCIS and CIS schemes with Ur = 2 re- 
lays. All techniques employ RLS algorithms for the estimation 
of the coefficients of the channel, the receive filters and the 
power allocation for each user (JPAIS only). The complexity 
of the proposed algorithms is quadratic with the filter length of 
the receivers and the number of relays n^, whereas the optimal 
MMSE schemes require cubic complexity. From the results, 
we can verify that the proposed adaptive joint estimation 
algorithms converge to approximately the same level of the 
MMSE schemes, which have full knowledge of the channel 
and the noise variance. Again, the proposed JPAIS-GPC is 
superior to the JPAIS-IPC but requires higher complexity and 
joint demodulation of signals, whereas the JPAIS-IPC lends 
itself to a distributed implementation. 

The next experiment considers the average BER perfor- 
mance against the normalized fading rate fdT (cycles/symbol). 



Fig. 6. BER performance versus number of symbols. The curves for the 
adaptive schemes are in solid lines, whereas those of the optimal MMSE 
schemes are in dotted lines. Parameters: \t = = 0.025 (for MMSE 
schemes), a = 0.998 (for adaptive schemes). 



N=16, SNR = 12 dB, K=6 users, P=1500 




Fig. 7. BER performance versus fdT. The parameters of the adaptive 
algorithms are optimized for each fdT. 



as depicted in Fig.|7] The idea is to illustrate a situation where 
the channel changes within a packet and the system transmits 
the power allocation vectors computed by the proposed JPAIS 
algorithms via a feedback channel. In this scenario, the JPAIS 
algorithms compute the parameters of the receiver and the 
power allocation vector, which is transmitted only once to the 
mobile users. This leads to a situation in which the power 
allocation becomes outdated. The results show that the gains of 
the proposed JPAIS algorithms decrease gradually as the fdT 
is increased to the BER level of the existing CIS algorithms for 
both rir = 2 and nr = 4 relays, indicating that the power al- 
location is no longer able to provide performance advantages. 
This problem requires the deployment of a frequent update 
of the power allocation via feedback channels. Therefore, 
the proposed JPAIS algorithms are suitable for scenarios with 
fdT < 10~^ for which there are performance gains. 

The algorithms are now assessed in terms of the mutual 
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information between the fc-th user and the base station as 
suggested in ifTTI and the normalized throughput (NT) defined 
as NT = R{1 - BER^ '"^^^^ in bits/time slot, where R 
is the normalized rate, M is the number of points in the 
constellation and P is the packet size in symbols. Since the 
protocols operate at full rate in a synchronous system and 
QPSK modulation is used the parameter used are i? = 1 and 
M = 4. The results are shown in Fig. |8] and indicate that 
the JPAIS algorithms can obtain gains in terms of the mutual 
information and the NT for a sufficiently high SNR level. 
When the SNR is low the use of more phases of transmission 
(or time slots) can degrade the NT and the mutual information 
of the system. The use of the proposed JPAIS algorithms 
with multihop transmission is not recommended in these 
situations. However, as the SNR is increased the proposed 
JPAIS algorithms obtain the best results with = 1,2. The 
JPAIS-GPC algorithm achieves the best NT results followed 
by the JPAIS -IPC and the CIS techniques. 



N=16, SNR = 12 dB, K=6 users, P=1500, f T=10"- 
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Fig. 9. BER performance versus f^T. The parameters of the adaptive 
algorithms are optimized for each /^T. 
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Fig. 8. (a) Mutual information (bits/Hz) versus SNR (dB) and (b) the 
normalized tliroughput (bits/time slot) versus SNR (dB). 



The last experiment, shown in Fig.|9] illustrates the averaged 
BER performance versus the percentage of errors in the 
feedback channel for an uplink scenario. Specifically, the 
feedback packet structure is employed and each coefficient 
is quantized with 4 bits. Each feedback packet is constructed 
with a sequence of binary data (Os and Is) and is transmitted 
over a binary symmetric channel (BSC) with an associated 
probability of error Pg. We then evaluate the BER of the 
proposed JPAIS and the existing algorithms against several 
values of the Pp. The results show that the proposed JPAIS 
algorithms obtain significant gains over the existing CIS 
algorithm for values of < 0.1%. As we increase the rate 
of feedback errors, the performance of the proposed JPAIS 
algorithms becomes worse than the CIS algorithms and are no 
longer suitable. This suggests the use of error-control coding 
techniques to keep the level of errors in the feedback channel 
below a certain value. 

VIII. Concluding remarks and extensions 

We have presented in this work joint iterative power allo- 
cation and interference mitigation techniques for DS-CDMA 



networks which employ multiple hops and the AF cooperation 
strategy. A joint constrained optimization framework and algo- 
rithms that consider the allocation of power levels across the 
relays subject to global and individual power constraints and 
the design of linear receivers for interference suppression have 
been proposed. A study of the proposed optimization problems 
has been carried out and has shown that the convexity of 
the problem can be enforced via an appropriate choice of 
the global and individual power constraints. A study of the 
requirements of the proposed and existing algorithms in terms 
of computational complexity and feedback channels has also 
been conducted. The results of simulations have shown that 
the proposed JPAIS techniques obtain significant gains in 
performance and capacity over existing non-cooperative and 
cooperative schemes. The proposed JPAIS algorithms can be 
employed in a variety of wireless communications systems 
with relays including multiple-antenna, orthogonal-frequency- 
division-multiplexing (OFDM) and ultra-wide band (UWB) 
systems. Prior work on asymptotic results has been reported in 
1221 and 1231 for 2-phase systems without power allocation. 
A possible extension would be to use the power allocation 
strategy of this work and investigate the asymptotic gains. 
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